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LONG-TERM  GOALS 

Our  long  term  goal  is  a  quantitative,  mechanistic  and  predictive  understanding  of  the  dynamics  of 
bubbles  and  bubble  populations  in  marine  sediments.  We  believe  that  this  information  can  be  used  to 
improve  and  test  acoustic  backscatter  models  for  sediments  and  to  better  understand  the  ebullitive  flux 
of  methane,  an  important  “greenhouse  gas”,  to  the  atmosphere. 

OBJECTIVES 

The  immediate  objective  is  to  develop  a  working  model  for  the  growth  of  a  single,  isolated  bubble  in  a 
marine  sediment.  This  modeling  is  being  guided  by  data  describing  bubble  growth  behavior  obtained 
in  laboratory  studies. 

APPROACH 

We  have  a  strategy  of  coordinated  laboratory  and  modeling  research  to  achieve  our  objective.  The 
laboratory  work  (directed  by  Bruce  Johnson  and  assisted  by  Bruce  Gardiner  and  Regine  Maass)  is 
focused  on  developing  an  understanding  of  the  growth  behavior  of  bubbles  in  sediments  through 
determining  the  relationship  between  stress  and  strain  during  bubble  expansion.  The  modeling 
research  (Bernard  Boudreau  and  Bruce  Gardiner),  guided  by  the  laboratory  results,  is  directed  toward 
developing  appropriate  model(s)  for  observed  bubble  growth  characteristics. 

WORK  COMPLETED 

During  the  third  (final)  year  of  this  granting  period,  we  have  made  substantial  progress  in  achieving 
our  objective.  We  have  characterized  the  mode  of  growth  of  bubbles  in  samples  of  marine  sediment 
as  growth  by  fracture  and  have  determined  values  of  the  important  variables,  including  the  critical 
stress  intensity  factor,  Kj^.,  and  Young’s  modulus.  These  results  have  been  described  in  a  paper 
submitted  to  the  Journal  of  Marine  Geology.  A  diffusion-reaction  model  that  incorporates  the 
experimental  results  for  the  mechanical  response  of  sediments  to  bubble  growth  has  been  developed  in 
which  a  distributed  source  of  methane  drives  bubble  growth  by  fracture.  This  work  has  been 
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described  in  two  papers,  one  submitted  to  Geochimica  et  Cosmochimica  Acta  and  the  other  to  the 
International  Journal  of  Heat  and  Mass  Transfer. 

More  recently,  we  have  taken  another  step  toward  achieving  our  long  term  goals  by  conducting 
experiments  into  the  dynamics  of  bubble  rise  in  sediments.  We  have  developed  the  methods  and  the 
apparatus,  and  are  presently  collecting  data  on  rise  rates  as  a  function  of  bubble  size  and  the 
properties  of  natural  and  surrogate  sediments. 

RESULTS 

A.  Experimental  Work,  Through  X-ray  images,  stress-strain  results  and  observations  with  surrogate 
materials  we  have  demonstrated  that  the  principal  mechanism  of  bubble  growth  in  our  sediment 
samples  is  by  fracture.  This  growth  mechanism  also  explains  coin  or  disk-shaped  bubbles  that  have 
been  reported  in  X-ray  tomography  images  in  other  natural  sediments  (Abegg  et  ah,  1994).  We  have 
shown  that  experimental  results  for  stress  and  strain  that  we  obtained  for  bubble  growth  in  natural 
sediment  samples  are  consistent  with  the  principles  of  linear  elastic  fracture  mechanics  (LEFM),  and 
we  have  used  these  experimental  results  to  determine  the  critical  stress  intensity  factor,  K^^.,  which  is 
the  principal  material  property  that  determines  fracture  behavior.  Values  of  for  our  sediment 
samples  (Cow  Bay,  Nova  Scotia)  range  from  about  2.3  XIO M  N  m'^^^  to  about  4.1  XIO  M  N 
In  addition,  we  have  estimated  the  critical  stress  intensity  factor  for  Eckernford  Bay  samples  by 
determining  the  eccentricity  of  bubbles  in  published  images  (Abegg,  et  al.,  1994).  The  obtained 
in  this  way  is  similar  to  our  Cow  Bay  results  and  is  about  5.5  XIO"^  M  N  m'^^.  The  manuscript 
describing  this  work  has  been  submitted  to  the  Journal  of  Marine  Geology. 


Figure  1:  Bubble  rising  by  buoyancy  driven  fracture.  Left  image  is  front  view  showing  bubble  5  cm 
long  by  3  cm  wide.  Right  image  is  side  view  showing  bubble  thickness  of  about  1  mm. 


In  studies  of  bubble  rise  in  natural  and  surrogate  sediment  materials,  we  have  determined  that  the 
mechanism  of  rise  is  by  buoyancy-driven  fracture.  Bubbles  rising  in  this  manner  typically  acquire  a 
shape  that  is  round  on  the  leading  edge,  more  square  on  the  trailing  edge  and  have  a  thickness  of  1  to 
2  mm.  We  are  now  conducting  experiments  to  determine  rise  speed  as  a  function  of  parameters  such 
as  Kij.,  E,  and  Ap.  We  then  plan  to  use  these  data  to  develop  a  model  of  rise  by  fracture  in  a  range  of 
sediment  types. 
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B.  Theory/Modeling.  We  have  used  the  experimental  results  deseribed  above  to  guide  development 
of  models  of  bubble  growth  by  fracture  with  methane  supplied  from  a  distributed  source.  Two 
models  have  been  developed  for  oblate  spheroidal  bubbles:  one  model  with  constant  bubble 
eccentricity  (DLCE,  i.e.,  diffusion  limited  constant  eccentricity)  and  the  other  with  eccentricity 
varying  according  to  linear  elastic  fracture  mechanics  (LEFM).  The  predictions  of  each  of  these 
models  have  been  compared  to  the  predicted  rate  of  growth  for  a  spherical  bubble.  Both  the  DLCE  and 
LEFM  results  show  that  oblate  spheroidal  bubbles  grow  substantially  faster  than  spherical  bubbles 
under  the  same  conditions  (figure  2)  and  can  reach  sizes  observed  in  marine  sediments  in  tens  of  days. 
In  additon,  the  LEFM  model  results  show  many  of  the  features  of  the  pressure-volume  curve  obtained 
for  injection  of  gas  into  natural  sediment  samples,  e.g.,  the  sawtooth  pattern  of  pressure  rise  and  fall  as 
gas  is  injected. 
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Figure  2:  Graph  showing  model  results  for  bubble  growth  rates  predicted  for:  oblate  spheroidal 
bubble  of  constant  eccentricity  (DLCE;  eccentricity  of  0.036  );  spherical  bubble  (DLCE; 
eccentricity  of  1 )  and  oblate  spheroidal  bubble  with  eccentricity  determined  by  the  principles  of 

linear  elastic  fracture  mechanics  (LEFM). 

[Both  oblate  spheroidal  bubbles  grow  substantially  faster  than  the  spherical  bubble;  the  LEFM 
bubble  grows  more  slowly  than  the  constant  eccentricity  bubble  to  about  5  mm  equivalent  radius 

and  faster  at  larger  sizes.] 


In  the  LEFM  treatment,  bubbles  at  small  sizes  grow  relatively  more  slowly  than  predicted  by  DLCE, 
but  at  larger  sizes  grow  faster.  One  particularly  interesting  prediction  of  the  LEFM  model,  and  a 
consequence  of  bubble  growth  in  an  elastic  medium,  is  that  small  bubbles  may  stop  growing  when  the 
bubble  pressure  required  for  fracture  exceeds  the  gas  pressure  in  the  pore  water.  Under  such 
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conditions,  larger  bubbles  grow  at  the  expense  of  smaller  bubbles.  The  manuseript  deseribing  this 
work  has  been  submitted  to  Geoehimiea  Cosmoehimiea  Aeta,  and  a  more  general  treatment  submitted 
to  the  International  Journal  of  Heat  and  Mass  Transfer. 

IMPACT/APPLICATIONS 

Bubbles  in  sediments  ean  seriously  eompromise  aeoustie  sensing  of  naval  mines,  destabilize 
struetures  that  rest  on  the  bottom,  and  transport  methane,  a  potent  greenhouse  gas,  to  the  atmosphere. 
Thus,  understanding  bubble  formation  (and  movement)  in  sediments  eonstitutes  an  important  praetieal 
and  seientifie  problem.  Our  findings  will  provide  information  that  eould  help  remove/filter  bubble- 
produeed  aeoustie  interferenee,  prediet  meehanieal  stability  of  sediments,  and  plaee  limits  on 
estimates  of  the  flux  of  methane  to  the  atmosphere. 

TRANSITIONS 

We  are  not  formally  eooperating  with  any  partieular  ONR  funded  projeet,  but  we  hope  to  integrate  our 
study  with  work  being  done  in  the  Bubble- Aeousties  DRI  in  the  future. 
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